STABILITY FOR SOLUTIONS OF WAVE EQUATIONS 
WITH C 1 1 COEFFICIENTS 



MIKKO SALO 

Abstract. We consider the stable dependence of solutions to wave 
equations on metrics in C ' class. The main result states that solu- 
tions depend uniformly continuously on the metric, when the Cauchy 
data is given in a range of Sobolev spaces. The proof is constructive and 
uses the wave packet approach to hyperbolic equations. 



1. Introduction 

We consider the wave equation in R< x R", 

{ (D 2 -A(x,D x ))u(t,x) =F(t,x), 

u\t=o = f, 
d t u\ t=0 = g. 

Here A(x, D x ) = a 13 (x)D Xi D x . is a uniformly elliptic operator, satisfying 
a %j = a ji^ anc L a l ^(x)d^j > c|£| 2 for £ € R n . We assume that the functions 
a 13 are in the space C 1 ' (R n ), with norm ||a|| c &-i,i = X^|a|<fcll^ aa lli°°- 

The question investigated in this article is the stable dependence of the 
solution u on the metric (a 13 ). Intuition for the problem can be obtained 
from the simplest possible case, namely the one-dimensional wave equation 
with constant sound speed c. Given / 6 L 2 (R), the equation 

(d 2 -c 2 d 2 x )u(t,x) = 0, 
u(0) = /, d t u(0) = 

has the solution u(t,x) = \[f{x — ct) + f(x + ct)]. Thus, the solution at 
time t = 1 is obtained by translating / by c units in the positive and nega- 
tive directions. Since translation of L 2 functions is a uniformly continuous 
operation, we see that u(l, ■ ) depends uniformly continuously in L 2 norm 
on the sound speed c. A stronger result may be obtained if the initial data 
is smoother: if / G H 1 , then ||/(- - c) - /( • -c')|| L 2 < II V/|| L 2 \c - c'|, and 
u(l, • ) depends Lipschitz continuously in L 2 on the sound speed. 

Our main result is the following theorem, which shows that uniformly 
continuous or Lipschitz dependence are valid also for general wave equa- 
tions with C ' metrics. In this introduction, we state the theorem only in 
the case where an initial velocity g is present. In Section 7 we will give 
the straightforward extensions to cases where driving terms F and initial 
positions / are present, and also where A(x, D x ) is replaced by a divergence 
form or Laplace-Beltrami operator. 



i 
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Theorem 1.1. Suppose M is a large constant such that 

IkicM <M, a%^'>M- 2 |£| 2 , (i) 

t€[-M,M]. (2) 

If —1 < a < 2, then for each g G if" there is a unique weak solution 
it G C([-M,M\]H a+1 )nC 1 ({-M,M\]H a ) of the problem 

r (Df-A^D^Ht,^ =o, 

| u\ t =o = 0, (3) 

[ d t u\ t =o = g- 

If — 1 < a < 2, and if A = (a* J ), -B = satisfy (pQ) and ua, ub are the 
corresponding solutions, then for any e > there is 5 > such that 

||itA(*) — UB(t)\\ Ha +i < e whenever \\A — B||co,i < 5. (4) 

Further, if -1 < a < 1 and g G H a+1 , then 

|[nA(t) - u B (t)\\ Ha +i < C\\A - B\\ c o,i WgWn^ (5) 

where C depends only on M and n. 

Here, by a weak solution we mean a function u G C([— M, M]; H a+1 ) n 
C 1 ([— M, M]; which solves the equation in the sense of distributions, 
and satisfies the initial conditions in the vector-valued sense. 

Existence and uniqueness of solutions to linear wave equations is of course 
classical and can be established via energy estimates under quite general con- 
ditions. We refer to [8] for a comprehensive account and further references. 
A constructive method, valid in the setting of our theorem, for proving ex- 
istence and uniqueness of weak solutions was introduced in [6] using a wave 
packet approach. 

Stability estimates such as the ones in Theorem II. II are also classical, see 
[S] and references given there. In Section 7 we give an easy argument which 
uses just the existence and uniqueness of solutions, without any precise 
knowledge about the solutions. 

The novelty is that our proof of Theorem 11.11 is constructive: an explicit 
expression for the solution is given, and the stability properties are deduced 
from that. The constructive method gives the same intuition to stability as 
in the n = 1 case, namely that stability for solutions should be the same as 
for translating functions. We use the wave packet approach introduced in [6j, 
where the initial data is decomposed into wave packets, and an approximate 
solution to the equation is obtained by translating wave packets along the 
Hamilton flow. The stability is exactly governed by this translation. 

The main motivation for this study comes from inverse problems in seismic 
imaging, where the analysis of solutions of wave equations has a key role. 
Many existing results, see [3] for a survey, assume a linearization about a 
smooth sound speed cq and use microlocal analysis and calculus of Fourier 
integral operators (FIOs). 

There has been recent interest, see [3], in the practically more realistic 
case where cq is not smooth, and in this case few results are known. One 
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reason is that there is no calculus of nonsmooth FIOs. However, solution op- 
erators for nonsmooth wave equations are understood quite well due to the 
wave packet approach of [6] , and a more precise analysis of these operators 
is expected to lead to new results. This work is an attempt in this direc- 
tion. The related work [3] discusses propagation of singularities (in terms of 
concentration of wave packets) and application of wave packets in numerical 
computations. We remark that wave packets are the same as curvelets [T], 
[2], which have been introduced in image processing as an efficient way of 
representing functions with singularities on smooth curves. 

The proof of Theorem 11.11 is based on constructing an explicit solution 
operator for the nonsmooth wave equation, following the method of [6]. The 
idea is to localize the initial data to dyadic frequency shells, and write each 
localized piece as a superposition of wave packets at the given frequency. 
The action of the wave group on a wave packet is well approximated by 
translating the packet along the Hamilton flow, and this gives an approxi- 
mate solution operator for frequency localized initial data. These are added 
up to obtain a parametrix for the full equation. 

Actually, to handle the nonsmooth coefficients, at frequency level 2 k one 
truncates the coefficients to frequencies less than 2 fc//2 and uses the Hamilton 
flow for the truncated metric. The error terms resulting from this will add 
up to a bounded operator, and this error can be iterated away by solving a 
Volterra equation. 

The precise construction of solution operator will be a combination of 
methods in [6] and [7J and it does not involve new ideas. However, in view 
of the stability result, we need to give the construction in detail to see how 
the operator depends on the metric. The main outline is the same as in [6], 
which used a discrete wave packet frame. Some computations are simplified 
if one uses instead a continuous wave packet representation (i.e. the FBI 
transform, see [S]). This was used in [9], [TU], [TT] for wave packets based on 
the Gaussian. We will follow [7J which used instead wave packets compactly 
supported in frequency, a property which will keep the dyadic frequency 
annuli separated. 

The crux of the stability proof is Lemma 16.44 which considers the sta- 
bility of translating along Hamilton flow. Lipschitz stability for translation 
involves a loss of one derivative, and the main point in the proof is that there 
is a smooth deformation of two Hamilton flows obtained by deforming the 
corresponding metrics. For the full result, one needs to check that the several 
corrections required to obtain an exact solution operator, most importantly 
the Volterra iteration, do not affect the stability given by translation. 

We first prove Lipschitz stability with a loss of one derivative, and the 
uniform continuity is an immediate consequence. Since there is a loss of 
one derivative arising from the flow, we can afford to lose one derivative in 
certain other estimates as well. In this way, we get stability in terms of the 
norm of the metric instead of C ' norm. 

The plan of the paper is as follows. Section 2 contains some basic facts 
about Hamilton flows, and Section 3 introduces the FBI transform. In 
Section 4 we outline the construction of the solution operator, and Section 5 
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contains the details. The stability result, Theorem ll.il is proved in Section 
6. Section 7 discusses variations of the stability result. 

Notation. We write Dt = idf and D Xj = \d Xj . The gradient with respect 
to x or £ is denoted by d x or d%, and D = id. Throughout the paper, M 
will be a large constant so that (pQ) and ((2|) are satisfied. We write a < b if 
a < Cb where C > depends only on M and the dimension n. Also, a ~ b 
means that a < b and b < a. We write L p = L p (R n ), similarly for the I? 
Sobolev spaces H a , the spaces C 0,1 and C > , and the Schwartz space 5? . 
The mixed norm spaces are denoted by L^H" = L p ([—M, M]\ H a ), similarly 
C{H% = Cj([-M,M};H<*). 

Acknowledgements. Research partly supported by the Academy of Fin- 
land. Part of this research was carried out during visits at MSRI and at the 
University of Washington, and I wish to express my gratitude to these insti- 
tutions for their hospitality and support. I would like to thank Hart Smith 
for generous advice, and Gunther Uhlmann for suggesting the problem. 

2. Hamilton flow 

We record for later use some elementary facts related to Hamilton flows. 
Let x{£.) De a smooth cutoff supported in the unit ball with \ = 1 f° r 
|£| < 1/2. We define smooth approximations of the coefficients a lJ by 

a l l{x)= X {2- k/2 D x yi{x). 

Then is supported in frequency in {|^| < 2 fc / 2 }, and satisfies derivative 
bounds 

ll^Xix- <2f max (°>i Q i- 2 ). (6) 

Consider the Hamilton flow related to Df — A^{x, D x ), where A^(x, D x ) = 
d£(x)D Xi D Xj . It will be useful to define this in terms of the functions 

p±(s,0 = ± X {2- k/2 D x )^Mx~J). 

The Hamilton equations are given by 

x{t) = dtp{x(t),£(t)), 

i(t) = -d x p{x(t),£{t)), 

where p = p^, x = x^, and £ = £ fe . Here h-> (d£p(x,£), -d x p(x,£)) is 

a smooth vector field in T*R n = R™ x (R^ \ {0}), so given an initial con- 
dition (x(0),£(0)) = (y, rj) G T*R n , the Hamilton equations have a solution 
(x(t),£(t)) depending smoothly on t, y, rj at least for small time. 

It is well known that the solution (x(t),£(t)) exists for all time. For, if 
it exists for t in an interval I = (—to, to), then ^ for t G I, and 
h(t) = |£(i)| 2 satisfies by the homogeneity of p 

h(t) = 2£(t) • i(t) = -2£(t) • = m(t)fc(t) 

where m(t) = — 2w(i) • d x p(x(t),uj(t)) and = Solving the 

equation gives h(t) = eh m ( s ) ds h(0) . But \m(t)\ < 2sup xenn \ (jj \ = i\d x p(x,uj)\, 
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and this argument shows that 

e -t svLp\d x p(x,w)\ufQ\i < < e *osup|dxp(a;,£j)||^/Q\| 

when tel. For x(t) we have \x(t)\ < sup xeR n|^| =1 |d^(x, w)| for t S I, so 
|x(t) — x(0)| < to^ 1//2 sup|(i^p(x, u)\ for tel. Thus stays in a 

compact subset of T*R n when t e I, and we may extend the solution past 
the endpoints of / and to all time. 

The following consequence will be used frequently: if (x(t),£(t)) satisfy 
the Hamilton equations and |£(0)| ~ A, then 

for \t\ < M. 

We write xt,s = Xt's f° r ^ ne canonical transformation (y, rj) i— > £(£)), 
where (x, £) solve the Hamilton equations with initial condition (x(s), = 
(y,T]). Then xt,s is a symplectic diffeomorphism of T*R n with inverse Xs,t- 

3. Wave packet representation 

Let g be a real, even Schwarz function in R n with ||<7||z,2 = (2ir)~2 , and 
assume g is supported in the unit ball. For A > 1 and y,x,£ € R n define 

g x (y;x,0 = ^e^<y^g(\h(y-x)). (7) 

This is a wave packet at frequency level A, centered in space at x and in 
frequency at £. The Fourier transform is given by 

g x ( V ;x,0 = X-^e-^g(X-H V - 0)- (8) 
The FBI transform of a function / E ^(R n ) is given by 



T\f(x,0= I f(y)gx(y;x,C)dy. 

Suppose A > 2 6 . Then if / is supported in |A < |£| < A, then T\f vanishes 
unless ±A < |f | < 2A. If F e J^(Rjf), the adjoint has the form 

T A F(y) = JJ F(x,09x(y,x,0dxd£. 

It follows that T A T A = I, and ||TA/|| i2{R 2n j = ||/|| L 2 (R n). 

The following result, stating the 1? boundedness of FBI transform type 
operators, is from [7J Lemma 3.1]. 

Lemma 3.1. Let g x £ be a family of Schwartz functions whose Schwartz 
seminorms are bounded uniformly in x and £. Then the operator T, defined 
for Schwartz functions by 



Tf(x,0 = j f(y)(g^)x(y;x,Ody 

is bounded from L 2 (R n ) to L 2 (R 2n ), and the adjoint T* given by 

T*F(y) = jj F(x,Z)(g Xt z)x(y\x,Z)dxd£ 

is bounded from L 2 (R 2n ) to L 2 (R n ). The norms of T and T* are bounded 
by a sum of finitely many Schwartz seminorms of the g x ^. 
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4. Outline of construction of solution operator 

We will outline the construction of a solution operator S(t) : g \— > u(t, ■ ) 
for the problem ([3|) . More details are given in Section [5j 

Let us start with the standard Littlewood-Paley frequency localization. 
Let x{0 be a smooth cutoff supported in the unit ball with x = 1 for 
|£| < 1/2, and take (5 k {D) to be a Littlewood-Paley partition of unity with 

oo 

/%(£) + Z>(£) = i 

k=l 

where (5q is supported in the unit ball, f3\ is supported in < |£| < 2}, and 
/?i(2- fc+1 6- 

We let g k = [3 k (D)g, and for each k we want to find an approximate 
solution of 

f (D?-A k (x,D x ))u k =0, 

| Ujt(O) = 0, 9iii fc (0) = g k . 

This will be done by reducing matters to first order hyperbolic equations. 
Consider the pseudodifferential symbol 

p±(x,0 = ±x(2- fc/2 ^)7I^) 

so that p^ is positively homogeneous of degree 1 in £ and smooth when 
£ 7^ 0, and satisfies for |£| = 1 

|d£afp£(x,OI <2f max (°'l a l- 2 ). 

Thus pfPk(0 £ ! / 2 - The Fourier transform in x, (pf) (•,£), is supported 
in {\r]\ < 2 k l 2 }. Since 

(p(x,D)/) A (77) = (2tt)-" / p(T/-C,0/(0de. 



we see that pf(x,D)f will be supported in | -yy ] ~ 2 fc if / is supported in 
1 77 1 ~ 2 fe . We will also use the symbols 

q± = X (2- k / 2 D x )(l/p±) 

which are homogeneous of degree — 1 in £, and q^f3 k (^) £ ^11/2- 
The approximate solution operators related to Q are 

E±(t)g = TtU±(t)T k (±Q±p k (D)g) 

where T k = T\ with A = 2 k , and U, (t) is translation for time —t along the 
Hamilton flow of p k , that is, 

U k ± (t)F = Fo X k Q f. 

We are ready to write down the first approximate solution operator for ([3]). 
It will be given by 

S(t)g = t ^ g k + J2( u t + U k~) 

fc<fco k>ko 
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where ko is a constant depending on M which will be chosen later, and 
where 

u t = E t(t)9- 

Since is of order —1, it is not hard to see that S(t) is an operator of 
order —1, i.e. it maps H a (R n ) to H aJrl (R n ). The main point is that the 
operator 

R±(t) = (D t + P±(y,D y ))E±(t), 
which is a half- wave operator applied to E, (t), is of order —1. This implies 
that when one applies the wave operator Df — A(x, D x ) to S(t), the resulting 
operator will be of order 0, which is one order lower than expected. Therefore 
we really have an approximate solution operator. 

However, since dtS(t)g\ t= o is only close to g but not equal to g in general, 
we need to correct the initial values of the operator. Thus instead of S(t) 
we use 

S(t) = S(t)(I + K)- 1 

where K is given by 

K = E (°) + ^fc (°) + \ R iPk(D) + l -R- k P k {D)) 

k>ko 

Here K^(3k(D) are order —1 operators given by 

If ko is chosen large enough then K will be small on H a (R, n ), and S(t) will 
be an operator of order —1 with S(t)g\t=o = and dtS(t)g\t=o = g- Thus 
S(t)g will be an approximate solution of Q with correct Cauchy data. 

It remains to show that one may iterate away the error and obtain an 
exact solution. To do this we seek a solution of ([3]) of the form 

u(t, x) = S(t)g(x) + / S(t, s)G(s,x) ds. 
Jo 

Here S(t, s) = S(t — s) is the operator corresponding to S(t) but where the 
initial surface is {t = s} instead of {t = 0}. Then S(t,s)g\t= s = and 
dtS(t, s)g\t= s = g, and one has 

d 4 2 (^ S{t,s)G(s,x)ds} =G(t,x) + J d?S(t, s)G{s,x) ds. 
We obtain 

{pi - A(x,D x ))u(t,x) = T(t)g(x) - G(t,x) + I T(t,s)G(s,x)ds 

Jo 

where T(t,s) = (D^ — A(x, D x ))S(t, s). As remarked above, we will show 
that T(t, s) is an operator of order 0, which in the present setting with a 
C 1 ' wave operator will mean that it is bounded on H a (R n ) for —1 < a < 
2. Then u will be a solution provided that G(t,x) = V(T(t)g(x)), where 
G = VF solves the Volterra equation 

G(t,x)- I T(t,s)G(s,x)ds = F(t,x). 
Jo 
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Since T(t, s) is bounded on H a also V is bounded on L^°H X , for — 1 < a < 2. 
Thus the full solution operator for (|3|) will be 

S(t)g(x) = S(t)g(x) + f S{t, s)V(T(s)g(x)) ds. 
Jo 

We will now give the details. 

5. Details of construction of solution operator 

Let (a*- 7 ) be a symmetric matrix of C 1, functions satisfying (pQ). We take 
M large enough so that ([1]) holds also for the truncated metrics (ou), and 
we also assume ([2]). We will use all the notations in Section [H Also, ko will 
be a sufficiently large integer, depending on M and n. 

We start by noting that if the Cauchy data is localized near frequency 
2 k , then the operators constructed in the preceding section preserve this 
localization. 

Lemma 5.1. Suppose the Fourier transform of / is supported in |£| ~ 2 k . 
Then T k f, P±f, Q±f vanish unless |f| ~ 2 fc , provided k > k . If F(x,£) 
vanishes unless |£| ~ 2 k , then T^F and U^(t)F vanish unless |£| ~ 2 k , 
provided k > ko. 

Using this lemma, it will be enough to consider a fixed frequency and the 
result will follow by summing over dyadic annuli. Thus, let A = 2 k , and 
write p = p^, U{t) = Uj^(t), etc. The main idea is that the wave evolution 
of a wave packet at frequency A can be well approximated by transport along 
the Hamilton flow. The transport operator will be L x ^ = 1st, given by 

Lx,£ = d$p(x, £) • D x - d x p(x, £) • D$. 

It is easy to see that U(t)F will solve the corresponding transport equation. 

Lemma 5.2. If F € ^(R™ x R^) then U{t)F satisfies 

(A + L x ^)U(t)F = 0, 
U(t)F\ t=0 = F 

Let now / be a function localized near frequency A. We write / as a 
superposition of wave packets at frequency A, 

f(y) = T{T x f{y) = ff T x f{x,£)g x {y;x,a)dxd(i. 

To get an approximate solution u to the half-wave equation (Dt+P y )u(t, y) = 
with tt(0, y) = f(y), where P y = p(y,D y ), we transport the wave packets 
for time t along the Hamilton flow of p. Then u is given by 

u(t,y) = J J T x f(x,^)gx(y,Xt,o(x,^))dxd^. 

Using that \t,o is a symplectic map with inverse xo,t, u will be equal to 

u{t,y) = T{U{t)T x f. 
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To measure how far u is from an exact solution, we compute 
D t u = T* x D t U(t)T x f = -T* x L x ^U{t)T x f 

= ff (U(t)T x f)(x,OL x ^gx(y;x,Odxd^ 
the last equality by integration by parts. Thus 

(D t + P y )u= JJ{U{t)T x f){x,0{Py + L Xii )g x {y-x,C)dxdi. (10) 

The following lemma, corresponding to Lemma 3.2], will be crucial. 
Lemma 5.3. Suppose |£| ~ A. Then 

(P y + L x ^)g x (y;x,0 = (g x> z)\(y;x,£) 

where g x £ is a family of functions whose Fourier transform is supported in a 
ball of radius 2 and the Schwartz seminorms of g x e are uniformly bounded 
in x, £. In fact g x ^(z) = m Xi ^(z, D z )g(z) where 

m x Az, C) = / (1 - s)d 2 s (p(x + s\- l ' 2 z, £ + sX^Q) ds. (11) 



Here Tn x ,£ satisfies symbol estimates uniform in x and £, 

\dy^m x ^z,C)\<(z) 2 , |C|<2. (12) 

Proof. We define g x ^ by the relation (g x ^)\ = (P y + L x ^)g\. Recalling the 
formula (J7|) for g\, it follows that 

(gx,ti)\ = (P y ~ idtp ■ d x + id x p ■ d^)gx = {P y + id^p ■ d y - d x p ■ (y - x))g\ 

= (2tt)-" j e*™(p(y, r?) - d^p(x, ■ r, - d xP (x, £) ■ (y - x))(g x y(ri) dr? 

= (2ir)- n J p(y,rj)-p(x,0-d^p(x,Oiv-0-dxP(x,Oiy-x) (9xT(v)dri 
since £; ■ d^p = p by homogeneity. By ([8]) , we have 

g x4 (z) = (2vr)-"| e iz < [pix+X^z, Z+\V\)- p ( x , £)-d x p{x, 6-A" 1/2 z 

-d^(x,0-x 1/2 c]g(0d<: 

which shows that {g x ,i) = outside a ball of radius 2. From Taylor's 
formula we see that the term in brackets is equal to m x ^(z, £) given by (jllH . 
If (x, £) = {x + sX~^ 2 z, £ + sA 1/,2 C), we compute 

9s(p(x,i)) = '^2[d XjXk p(x,i)X~ 1 z j z k + d x .£ k p(x,£)zj( k + 9^ k p(x,i)XCj(k] 

(13) 

which gives (I12p . □ 



Now we can prove that the half-wave operator applied to the approximate 
solution operator gives an operator of order 0. 
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Lemma 5.4. Suppose / is supported in |£| ~ A. Then 
\\(D t + P y )T* x U(t)T x f\\ L i<\\f\\ L 2. 

Proof. If / is supported in |£| ~ A, the same holds for U(t)T\f. The result 
follows from (|1U[) , Lemma 15.31 and Lemma 13.11 □ 

Next we apply a second order wave operator D 2 — P 2 to the approximate 
solution operator, and show that there is a loss of one derivative. 



Lemma 5.5. Suppose / is supported in |£| ~ A. Then 
\\(D 2 -P 2 )T* x U(t)T x f\\ L2 <\\\f\\ L i. 
Proof. We write u = T*U(t)T x f and 

(A 2 " Py)u = (D t + Pyfu - 2P y (D t + P y )u. 
Since u and (Dt + P y )u are supported in |£| ~ A and p is of order 1, we have 

\\P y (D t + P y )u\\ L2 <X\\f\\ L2 

by Lemma 15.41 

Using ((HIl) and writing (P y + L x ^)g x = (g x ^)x, we get 

(D t + P y ) 2 u = ff (U(t)T x f)(x,0(Py + L x ^(g x ^ x (y;x,Odxd^. 

Since 

D x (g x ^)\ = -D y (g x ^)x + (D x 9x^)x, 

D^{g x ,^)x = (y - x)(g x ^) x + (D^g x ^) x , 

the argument in Lemma 15.31 gives 

(P y + L x ^)(g x ^) x = (m x> z(z, D z )g x ^) x + (L x ^m x ^(z, D z )g) x . 

The Schwartz seminorms of g x £ are uniformly bounded in x, £, and (|12|) 
shows that the same applies to the Schwartz seminorms of m x ^(z, D z )g x ^. 
Lemma 13.11 gives 

(U(t)T x f)(x, £)(m x £(z, D z )g X £)x{y;x,£)dxd£\\ L 2 (:R n ) < \\f\\ L a 

It remains to study the symbol seminorms of m x ^(z,() = L x ^m x ^(z, £). 
From (1111) we obtain 



rh x ^z, C) = f (1 - s)d 2 Up(x, • D xP (x + sX'^z, i + s\ l ' 2 
Jo L 

- d x p(x, ■ D^p(x + s\~ 1/2 z, £ + sA 1/2 C)j ds. (14) 
Since |£| ~ A a computation as in (fT3|) shows 

\^m X;i (z,C)\<\ 1/2 (z) 2 , |C|<2, (15) 
and Lemma 13 . 1 1 implies 

|] y I* (t/(i)T A /)( a ; ) e)(m^(^, J D.)5)A(3/;^,Od a; ^|| i2(R .™ ) < aV^H/H^. 
This concludes the proof. □ 
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We move on to the full parametrix, where all the frequencies are added 
up. The errors resulting from truncation are handled as in [6, Theorem 4.5]. 

Lemma 5.6. S(t) is a bounded operator H a — > H a+1 , for all a £ R. If 
-1 < a < 2 then (Df - A(x,D x ))S(t) is a bounded operator H a -> H a . 

Proof. Since qj: is order —1 one has E,(t) : H a — > H a+1 , and since E^{t)g 
is localized near |£| ~ 2 k the sum converges in H a+1 and : H a — ► H a+1 . 
Let now — 1 < a < 2. We write 

(A 2 " ^ D x ))S(t)g = ~tJ2 A (x, D x )p k (D)g 
+ (A 2 " D x ))(u+ + u^) + £ (M%, D x ) - A(x, D x ))(u+ + u^) 

k>ko k>ko 

and write the last expression as T\(f)g + T2(t)g + T%{t)g. It is clear that 
\\Ti(t)g\\ H <* < Hffllif"- For T 2 (t) we write 

Z) t 2 - A k (x, D x ) = D 2 t - (P±) 2 + (P±) 2 - A k (x, D x ). 

One has 

||^(A 2 -(p+) 2 KIIh«<NI^ 

k>ko 

by Lemma 15.51 From (jGj) we obtain 

(^-fMKle^ (16) 

which shows that (-P^) 2 — ^4fc(x,-D x ) is of order 1 and 

\\J2((Pk) 2 -Mx,D x ))u+\\ Ha < \\g\\ H a. 

k>ko 

Similar results hold for u^. 

For T^t) we need to show that 

r : i-> ^ (afc(x) - a(x))D Xi D Xj u+ 

k>ko 

is bounded on where a& = and a = a 13 . We will write (3j{D)T = 
J2k>k r jkD XlXj u+ where 

r iJfc = /9 i (£))(o(a:)-o fc (x))^ fc (£>). 

Here $k(0 are cutoffs to |£| ~ 2 fc which are 1 on the frequency support of 
u+. We let l = l (M) be an integer so that supp(/3 fc (£)) Q {2~ lo+k < \£\ < 

By looking at the supports on the Fourier side, we get 

Pj(D)(a - a 2k -4i )Pk(D), j <k- 3l , 
>=<( pj(D)(a-a k )p h (D), k - 3/ < j < k + 3Z , 

0j(D){a - a 2j _ 4lo )p k (D), j>k + 3Z . 
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Using that \\a — a^H^oo ^ 2~ fc we obtain 

E ^ 2k \\ 9k \\h+ E 2-^2^11^11^ 

k<j-3l \j~k\<3l 

+ £ 2- 4fc 2 2fc || 5fc ||| 2 . 
fe>i+3/ 



It follows that r is bounded on H a with \a\ < 1. We compute 

iXjX[ 



D Xl Tg = ^2 {Dia k )D XrX] ul - D L a ^ D XiXj u^ + E ~ ") l) >-. r ■ 



k>ko k>ko k>ko 

and note that this is bounded H a+l — > i7 a if \a\ < 1, due to the C > 
regularity of a and the argument above. This shows boundedness of T for 
-1 < a < 2. □ 

In the next lemma we correct the value of dtS(t) at t = 0. 
Lemma 5.7. If &o is sufficiently large, then 

K = E (°) + ^* (°) + + \RkfaiP)) 

k>ko 

has norm < 1/2 on H a . The operator 

S(t) = S(t)(I + K)- 1 
will satisfy 5(t)g|t=o = and d t S(t)g\ t= o = g- 

Proof. One has \\Rf(0)g\\ L 2 < 2- k \\g\\ L 2 by LemmaEl Also, R^(3k{D) is of 
order —1 by looking at the symbol expansion of (Pj^Q^ — I)Pk(D) and using 



(|16fl . The norms do not depend on any previous value of ho, and we may 
choose fco so large that K will have norm < 1/2 on H a for a £ [— aO) a o]> 
for any ao > 0. Thus I + K will be invertible on these spaces and the norm 
of the inverse will be < 2. 

Since S(0)g = 0, also S(0)g = 0. The derivative is 

d t S{t)g = 9k + i E ( D t u t + A«fc )• 

k<ko k>ko 

We write D t uf = R^(t)g - P^u^. Since u^(0) = ^Q^gk, we get 
^5(t) 9 | 4=0 = E 5 fc + E (»^fc (% + (% + ^ + 

fc<feo k>k(t 

It follows that d t £(%| t =o = (I + □ 

It now remains to show that one obtains a full solution operator from S(t) 
by introducing a correction by solving a Volterra equation. The solution 
operator for the Volterra equation is the following. 

Lemma 5.8. Suppose T(t,s) is bounded on H a , with norm uniformly 
bounded by C when t, s € [-M, M\. There is a bounded map V on L^°H° 



r, 
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with norm bounded by e , such that for any F(t,x) € L^°H X , G = VF 
solves the equation 

G(t,x)-[ T(t,s)G(s,x)ds = F(t,x). 
Jo 

Proof. Define (VF)(t,x) by 

00 ft rsi fSj-i 

F(t,x)+} j / /•••/ T(t,si)T(s 1 ,s 2 ) ■ ■ ■T(sj-i,Sj)F(sj,x)ds j ■ ■ ■ dsi. 
■ =1 Jo Jo Jo 

It is easy to check that the series converges in LfH x to a solution which 
satisfies the desired norm estimate. □ 

As in the outline, we write S(t, s) = S(t—s) as the operator corresponding 
to S(t) but with initial surface {t = s}. If s, t G [— M, M] then S(t, s) has the 
same properties as S(t), except that S(t, s)g\t= s = and dtS(t, s)g\t= s = g. 
We let 

T(t, s) = {pi - A(x, D x ))S(t, s), T(t) = T(t, 0), 

so T(t,s) and T(t) are bounded on H a when — 1 < a < 2 by Lemma 15.61 
We let V be the solution operator for the Volterra equation corresponding 
to T(t,s). 

The argument in the outline shows that if 

S(t)g(x) = S(t)g(x) + /* S(t, s)V(T(s)g(x)) ds, 
Jo 

then {D 2 t - A(x,D x ))S(t)g = and S(t)g\ t=0 = 0, d t S(t)g\ t=0 = g. One 
only needs to check that the time derivatives are justified, but this may be 
done as in [6] and will not be needed for stability considerations. This ends 
the construction of the solution operator. 

6. Stability 

We now proceed to prove the stability part of the result. Let A{x) = 
(a l i(x)) and B(x) = {W{x)) be two symmetric matrices satisfying (jTJ) , and 
take M so large that the truncated metrics [cff) and (bV) also satisfy ([I]). 
Also assume that t satisfies ([2]). 

We write Pa = P^ for the operator at frequency A = 2 k defined in terms 
of the metric A, and similarly Qa, Ra, -Ra(*) etc. The following operators 
depend norm continuously on the metric. 

Lemma 6.1. If / is supported in |£| ~ A, then 

\\(P A -P B )f\\ L2 <\\\A-B\\ L ~\\f\\ L2 , 
\\(.QA-QB)f\\iP<X- 1 \\A-B\\ L c \\f\\ I ^, 
\\(RA-RB)f\\^<^ 1 \\A-B\\ L ^\\f\\ L2 . 

Proof. Consider h(x,£) = (A"/ 2 X (A 1/2 • ) * [F(A k ( ■ , £)) -F(B k ( ■ , £))])&(£), 
where f3\ is a cutoff to |£| ~ A and F(t) = t 1 / 2 . We wish to show that 

\d^h(x,0\ < C M ,ad A - B \\L°°^~ m+1 * H - (17) 
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This will show the estimate for Pa — Pb, and the estimate for Qa — Qb 
follows from a similar result with F(t) = i -1 / 2 . 

we let the x-derivatives hit the mollifier, which gives the desired 

growth. We may thus assume that a = 0. Each ^-derivative hitting 
gives A -1 , so we only need to look at the case when the ^-derivatives hit 
F{A k ) - F(B k ). We write 

F(A k ) - F(B k ) = f F'(rA k + (1 - r)B k ){A k - B k ) dr. 
Jo 

The matrix rA k + (1 — r)B k satisfies ([T]). Consequently 

\&l [(F(A k ) - F(B k ))p k (0]\ < CmA A ~ BWl^X 1 ^ 

as desired. 

If P, Q are pseudodifferential operators the symbol of PQ is 
„ dfpD^q „ f 

a{PQ) = Y, + E / e ~^M^ £, y, V) dy d V 

\a\<N ' \a+/3\=2N 

where the last terms are oscillatory integrals, and 

7<a,/3 J ° 

Note that \a + /3 - 7| > N. We have 

(R A - Rb)Px(D) = Pa(3\(D)QaI3\(D) - P B Px(D)Q B x (D) 
for a suitable (3\. Suppressing the cutoffs, this has the symbol 

/ £j r> \ V s df P AD* q A-d?PBD« qB 

a(R A - Rb) = PAQA ~ PbQb + > — 

z — ' od 

0<\a\<N 

+ E e ~ iy ri (. s £f3( x i^yiv) - s%p(x,£,y,v))dydri 

\a+P\=2N 

The principal symbol is paQA-PbQb = (1 ~x( 2 ~ fc/2 Ar))(V.PB - since 
q = x(2 _fc / 2 D x )(l/p) , and the arguments above show that the corresponding 
operator has norm < A 1 1 1 ^4 — -B||loo on 1? . It is easy to see that the terms 
with < \a\ < N have the same bound. Finally, if I7I + |<5| < 2n + 1 then 

[d^'af {si P - s^)(x, e, y, r,)\ < ca-t +n+i-\P'HW\ \\ A -B\\ L ~. 

Taking N large enough and using standard estimates for oscillatory integrals 
gives the L 2 bound for Ra — Rb- D 

Let gg^ be the Schwartz functions in Lemma 15.31 defined in terms of the 
metric A. 

Lemma 6.2. The Schwartz norms of g x ^—Qx^ are ~ A||^4 — -B||z/°° > uniformly 
in x and 
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Proof. It is enough to show that 

\d^(m^ - m^)(z, C)| < X\\A - B\\ L ~{z) 2 , \(\ < 2. (18) 

This follows from the expression (jlip for m X) £, the computation (|13p . and 
the estimate (fT7|) . □ 

Remark. It is also true that the Schwartz seminorms of cQ^ — are 
< \\A — B\\qi,i, which follows by replacing (|17p with the alternate estimate 

\d%dPh{x,£)\ < C M>a>0 \\A - BWc-^x 1 -^^ ^ 2 ) . 

However, due to translation along Hamilton flow, there is a loss of one 
derivative in the stability estimate in any case. Therefore we can afford to 
lose one derivative in other estimates as well. This results in stability in 
terms of the C ' 1 norm of the metric instead of C 1, . 

Lemma 6.3. \\(R A (0) - R B (0))f\\v> < U ~ B\\ L <»\\f\\i?. 
Proof. One has 

(R A (0)-R B (0))f =l( JJ (Tx(Q A -QB)Px(D)f)(x,0(g£s)x(y;x,Odxdt 
+ JJ (TxQ B Px(D)f)(x,0(g^ - g^)x(y;x,Odxd^j . 

The result follows from Lemmas 13.11 15.31 16.11 and 16.21 □ 

From the preceding results, and from the factorization 

(/ - KaY 1 -(I- Kb)- 1 = (I- K A )- l (K B - K A ){I - K B )-\ 

we see that (I+K)' 1 , the operator which corrects the initial values, depends 
norm continuously on the metric: 

+ Ka)' 1 - (I + KB^Wn^Ha < \\A-B\\ L oo. (19) 

Next we consider the stability of the flow operator. Here we will lose one 
derivative to get Lipschitz stability. 

Lemma 6.4. If / vanishes unless |£| ~ A, then 

\\n(U A - U B )Txf\\ L 2 < X\\A - Blloo,! ||/|| L2 . 

Proof. Fix t and write xa = Xot> Xb = Xot- One has 

Tl(U A - U B )T x f(y) = Jf [Txf( XA (x,0)-Txf(XB(x,0)]gx(y,x,OdxdC. 

We let C r (x,0 = rA x (x,0 + (1 - r)B x (x,£), and let $ r = X c r = (x r ,Cr) 
be the flow corresponding to the metric C r . Then ($r)r6[0,il i s a smooth 
family of symplectic diffeomorphisms of T*R™, and we have 

Txf(XA(x,0)-T x f(XB(x,0) = [' (d x ^T x f)(Mx, 0) ■ d r Mx,0 dr. 

Jo 
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Let h(s) = h(s,r,x,£) = (x r (s),£ r (s)/A) where (x r (0),£ r (0)/A) = (x,£/A). 
Differentiating the Hamilton equations for (x r ,£ r /A) with respect to r, and 
using 

a ' pc ' = x(x ' D ' ] W ,0 ' 

we obtain 

\(d r h)\s)\ < \d r h(s)\ + \\A-B\\ c o,i. 

Since d r h(0) = (0,0), Gronwall's inequality shows \d r h(s)\ < \\A — B\\(jo,i 
for |s| < M. This implies 

l^r| < 11-4 - B||cto,i, 
l^rCrl < \\\A-B\\ CO ,i. 

Using that d x .T\f = T x (d 3 f) and %T A / = X^ 2 f jjX f where f jA f(x,$ = 

(f, (izjg)\(-;x,0)i we s et 
T A *(f/ A - U B )T x f(y) = J J J [(9 r x r ),(x,0T A (5 J /)(cI> r (x,0)9A(y;x,e) 

+ (dr^(x,0^' 1/2 T j ,xf(M x ,0)9x(y;oc,o] dxd^dr. (20) 

Lemmas 13.11 and 15.31 imply the desired estimate. □ 

Lemma 6.5. Suppose / is supported in |£| ~ A. Write M A f = (D t + 
P A )T*U A T x f and N A f = (D 2 - (P A ) 2 )T*U A T x f . Then 

\\(M A - M B )f\\ L2 < X\\A - BIIcm ||/|| L2 , (21) 

\\(N A - N B )f\\ L2 < X 2 \\A - B\\ CO a \\f\\ L 2. (22) 

Proof. From (|10p we have 



(M A -M B )f = J]{U A - u B )T x f( x ,0(g£,t)(y,x,Odxdt 

+ 11 U B T x f(x,0(g^-g^)x(y;x,Odxd^ 



where g A ^ is as in Lemma 15.31 That lemma, the argument in Lemma IS~ 
and Lemma 16.21 give (|21j) . For the other estimate we write 

D 2 - P 2 = (A + Py)(D t + Py) ~ 2Py(D t + Py) 

which gives 

(N A - N B )f = [(D t + P A )M A - (D t + P y B )M B ]f 

- 2(P y A - Py)M A f - 2P y B (M A - M B )f. 

The last two terms have L 2 norms < A 2 ||^4 — ^Hc - 1 II/IIl 2 by Lemma 16.11 
and (|2ip . The first term has the form 

[(Dt+P A )M A -(D t +Py B )M B ]f =ff (U A -U B )T x f(x,i)(g A ^)( y] x^)dxdi 

+ J J U B T x f(x^)(g A ti -~g B tii ) x (y,x,i)dxdi 
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where g^ = m^(z,D z )g^ + rh^(z, D z )g, using the notation in Lemma 

15.51 The Schwartz seminorms of g^ are < A 1 / 2 by (fl"5l) . and those of 

9x,t ~ 9x£ are < x2 W A - B h°° by ([THD, CLU), (HZD- Now (|22|) follows in the 
same way as (|2ip . □ 

We proceed to prove stability results for the operators where all the fre- 
quencies are summed up. 

Lemma 6.6. \\(S A (t) - S B {t))g\\ H <* < \\A - B\\ c o,i\\g\\ H « ■ 

Proof. Because of (fT9l) it is enough to prove the estimate for S A (t) — S B (t). 
But we have 

(S A (t) - S B (t))g = J2 ((^U*) " K,B(t))g + (E~ A (t) - E^ B (t))g) 

k>ko 

where 

K A (t) - E± B (t))g = mU± A (t) - U± B (t)]T k C-Q± A (3 k (D)g) 

+ nut B {t)T k {^[Q% A - Q% B ]fa{D)g). 

By Lemmas 16. II and 16.41 one gets 

\\(E± A (t) - E± B (t))g\\ L 2 < \\A-B\M\g\\ L2 . (23) 
We sum up these estimates and use frequency localization to obtain 
\\(S A (t) - S B (t))g\\ Ha < \\A-B\\ c0 ,4g\\ Ha . 

□ 

The next lemma considers the error T A (t) = (D 2 — A(x, D x ))S A (t). 
Lemma 6.7. \\(T A (t) - T B (t))g\\ H °< < \\A - B\\ c o,i\\g\\ Ha +i when \a\ < 1. 
Proof. As above, it is enough to consider the operator 

f A (t)g = (A 2 - A(x,D x ))S A (t)g = -t J] A(x,D x )g k 

k<ko 

+ (A 2 - A k (x, D X ))(E+ A (t)g + E~ A {t)g) 

k>ko 

+ (M*> D *) - M*, D x )){El A {t)g + E- A {t)g). 

k>ko 

We write the last three terms as T A j(t)g for j = 1, 2, 3. Clearly 

\\(T AA {t)-f BA (t))g\\ Ha < \\A-B\\ c oa\\g\\ H «- 
For the second term we use that 

(A 2 - A k (x, D x ))E± A (t)g = [D 2 t - {P^ A f)E% A {t)g + C% A E±Jt)g 
where C^ A = (P^ A ) 2 - A k (x, D x ). The first term is just N A (^Qf A g k ), and 

II^Qt^fc) - N B { l -Q% B g k )\\ L , < 2 k \\A - B\\a>,i \\g\\v 
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by (|22p and Lemma 16. 11 Writing 

C%a ~ C% B = {Pt A ? ~ (Kb) 2 ~ (Mx,D x ) - B k (x,D x )), 
and using the argument in the end of Lemma 16.11 we see that 

IK^M " C k,B)h{D)g\\ L * < 2 2k \\A - B\\ L ~\\g\\v 
when $k(0 is a cutoff to |£| ~ 2 k . Lemma IB31 shows 

\\(f At2 (t)-f Bt2 (t))g\\ Ha < \\A-B\\ c o,i\\g\\ Ha+ i. 
For the last term we write 

(f Ai 3(t)-f Bj3 (t))5 = EE(^( x '^)-^(^^))(( s M(*)-^(*))^ 

± k>k 

± k>k 

The discussion in Lemma 15.61 and (|23|) give the bound for the first term. 
The second term is handled as in Lemma l5.6( except that we use 

\\a - a k (x) - (b- b k (x))\\ L oo < \\A - B\\ c o,i2- k / 2 . 

One first gets the bound for — 1 < a < and then for \a\ < 1 by computing 
the derivative. The result follows. □ 

We note that the last estimate holds for the operators T(t, s) with uni- 
form constants when t, s £ [—M,M]. The final estimate we need is for the 
Volterra solution operator. 

Lemma 6.8. \\(V A — V B )F \\l%°H!* < \\A — B\\ c o,i \\F\\ LOOH <x+i when \a\ < 1. 

Proof. Recalling the definition of V from Lemma 15.81 we need to estimate 
the Lf^H" norm of terms of the form 

ft fSl fSj-1 

I(t,x)= /•••/ T B (t,si)---T B (si^ 1 ,si)(T A -T B )(si,si +1 ) 
Jo Jo Jo 

T A (si +1 , s l+2 ) ■ ■ ■ Ta{sj-i, Sj)F(sj,x) dsj ■ ■ ■ dsi 

Choose C = C(M) such that for t, s £ [-M, M], 

\\T A (t,s)g\\ H c> + \\T B (t,s)g\\ H ° < C\\g\\ H «, -1 < a < 2, 
\\(T A - T B )(t, s)g\\ H a < C\\A - B\\c°a \\g\\ H ^ , -1 < a < 1. 

Using the H a estimate for each T B , the H a+1 — » H a estimate for T A — T B , 
and then the H a+1 estimate for each T A gives 

\\I{t, -)||jfa < -—\\A-B\\ c o,i\\F\\ Lrm+ i. 

3- * 

There are j terms of the form I(t,x) at level j. It follows that 

\\(V A - V B )F\\ L?HS < ( f; \\A - B\\ c0 ,i \\F\\ LrHS+1 . 

□ 
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For the full solution operator we have 

(S A - S B )(t)g = (S A - S B )(t)g + f\s A - S B )(t, s)V A (T A (s)g(x)) ds 

Jo 

+ [ S B (t,s)(V A -V B )(T A (s)g(x))ds 
Jo 

+ I S B (t,s)V B ((T A -T B )(s)g(x))ds. 
Jo 

Consequently, for g G H a with — 1 < a < 1, 

\\{S A ~ S B )(t)g\\ Ha +i < \\A - B\\ c oa \\g\\ H <*+i 

by Lemmas 16. 6| 16.71 and 16.81 This shows (0) • 

To prove ([!]), take g G H a with — 1 < a < 2, and write g = g s + g r with 
g s G if 2 and Hgyllfl"" small. Using the triangle inequality and (0), we get 

\\(S A (t) - S B {t))g s \\ H 3 < \\g s \\ H 2, 

\\(S A (t) - S B (t))g s \\ H 2 < \\A - B\\ c o,i\\g s \\ H 2, 

and by interpolation \\(S A (t) - S B (t))g s \\ H 3- K < \\A - B\\£, 0A \\g s \\ H 2. Thus, 
given e > 0, by choosing ||g r ||^ a and 5 > small enough we obtain 

\\(S A (t) - S B (t))g\\ Ha+1 < \\(S A (t) - S B (t))g s \\ Ha+ i 

+ \\S A {t)g r \\ Ha +i + \\S B (t)g r \\ Ha +i < e 
whenever ||^4 — SH^o,! < S. This ends the proof of Theorem ll.il 

7. Extensions 

Here we give some extensions of Theorem 11.11 following Section 4 of [6]. 
The solution operators in each case are obtained from S(t) similarly as in 
[6] , and the stability results follow from the arguments given in Section [6J 
Therefore we will omit the proofs in this section. We emphasize that for 
each theorem below, the stability proof is constructive. 

Theorem 7.1. Assume (P) - ([!]), and let —1 < a < 2. Suppose that 
/ G H a+1 , g G H a , and F G L\H^. Then there is a unique solution in 
C t °#" +1 n C\H% for the problem 

' (D? -A(x,D x ))u(t,x) =F(t,x), 

u\ t =Q = f, 

d t u\ t =o = g. 

The solution satisfies 

IK*)llir«+i < \\f\\ H -+i + \\g\\H° + \\ f \\l\hs- (24) 

Also, let A = (a lJ ), B = satisfy ([T]), and let u A , u B be the corresponding 
solutions. If —1 < a < 2, then for any e > there is S > such that 

— u B{t)\\H a + 1 < e whenever \\A — B\\ c o,i < 5. (25) 
Further, if -1 < a < 1 and / G H a+2 ,g G H a+1 ,F G L 1 iJ£ +1 3 then 
11^(0 -«b(*)II^+ 1 < P- 5 llc^ 1 (ll/ll^+ 2 + lbllH"+ 1 + ( 26 ) 
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Remark. Here we show how the existence and uniqueness part of Theorem 
17.11 can be used to prove stability of the map A i— > u(t, ■ ) in Theorem 11.11 
It is enough to prove ([5]) for — 1 < a < 1, and the uniform continuity will 
follow as in Theorem 11.11 Let A, B be two metrics satisfying ([1]), and let 
g G H a . We denote by ua, ub the solutions to ([3]) with data g, and write 
v = ua — ub ■ Then v satisfies 

f (A 2 - A{x, D x ))v = (A(x, D x ) - B(x, D x ))u B , 

{ v(0) = d t v(0) = 0, 

which implies 

\\v(t)\\ Ha +i < \\A - B\\ c o,i( sup 11^(^)11^+2) < \\A - B\\ c o,i\\g\\ Ha +i. 

\t\<M 

This is the required bound. 

Now we consider the wave equation with divergence form operator 

A D (x,D x )u = D Xi (a ij (x)D x .u). 

Theorem 7.2. Assume (pQ) - ([2]), and let -2 < a < 1. Suppose that 
/ £ H a+1 , g G H a , and F G L\H X . Then there is a unique solution in 
C^H^ +l n C}Hg for the problem 

f (D 2 t -A D (x,D x ))u(t,x) =F(t,x), 

u\t=o = f, 

d t u\ t =o = g. 



The solution satisfies (1211) . Also, solutions satisfy (|25l) and (|26|) . if the ranges 
for a are replaced by —2 < a < 1 and —2 < a < 0, respectively. 

Finally, we consider the Laplace-Beltrami type operator 

A L (x,D x )u = p{x)~ l D Xi {p{x)a ij {x)D Xj u) 

where p = (det (a 1 - 7 )) 1 / 2 . 

Theorem 7.3. Assume ([T]) - ©, and let — 1 < a < 1. Suppose that 
/ £ H a+ , g G iJ a , and i* 1 G L\H X . Then there is a unique solution in 
C t °tf£ +1 n C/if" for the problem 

(Df - A L (x, D x ))u(t, x) =F(t,x), 

u\t=o = f, 

d t u\ t =o = g. 

The solution satisfies (|24p . Also, solutions satisfy (|25p and (|26p . if the ranges 
for a are replaced by — 1 < a < 1 and — 1 < a < 0, respectively. 
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